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universe evolves to a state independently of the initial conditions at the edge 
of chaos. The critical state is an attractor of the dynamics. Random metric 
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fluctuations exhibit noise without any characteristic length scales, and the 
qv | power spectrum for the fluctuations has a self-similar fractal behavior. In 



the early universe, the metric fluctuations smear out the local light cones 



£ry removing the horizon problem. 
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One of the most important problems in modern cosmology concerns the fine-tuning 
necessary in the standard cosmology based on general relativity (GR). Why is the universe 
so close to being spatially flat after evolving for more than 10 gyr? Why is it so isotropic 
and homogeneous? How could such a critical state of the universe come about without a 
severe fine tuning of the parameters? 

The standard explanation for these questions has been the inflationary models 0. These 
models have faced problems that arise mainly from the need to fine tune certain parameters 
and initial conditions, e.g., the degree of inhomogeneity of the initial universe, or in Linde's 
"chaotic" inflation the need to fine tune parameters at the Planck energy. In the following, 
we shall study a self-organized universe which naturally evolves to a critical state without 
detailed specification of the initial conditions. The critical state is an attractor of the system 
which does not need to be fine tuned. This is in contrast to the inflationary models which 
have an attractor mechanism with critical phase transition points that need to be fine tuned. 
In statistical mechanics both kinds of attractor mechanism are known to occur for physical 
systems. In contrast to the inflationary paradigm, we shall be concerned with a self-organized 
universe in which the spacetime metric fluctuations [fj| have a high degree of cooperative 
effects at the edge of chaos. Recently, the mixmaster universe has been shown to be chaotic, 
i.e., it does not evolve as a self-organized system in the early universe |§. 

In a recent new approach to gravitational theory ||, it has been proposed that at some 
length scale much larger than the Planck length, l p « 10 -33 cm, the spacetime geometry 
is fluctuating randomly. In classical GR it is assumed that the spacetime manifold is C 2 
smooth down to zero length scales. This seems to be an unacceptably strong hypothesis 
considering that known physical systems possess dynamical noise at some length scale and 
that cooperative effects are know to occur for many systems proportional to V° and not 
V~ a (a > 0) where V is a characteristic volume for the physical system. The metric ten- 
sor was treated as a stochastic variable and for a given three-geometry 3 Q, a stochastic 
differential equation for the momentum conjugate variable was obtained. A Fokker-Planck 
equation was derived for the probability density leading to statistical mechanical predictions 



for gravitational systems. 

Spatially extended dynamical systems with both temporal and spatial degrees of free- 
dom are common in biology, physics and chemistry. These systems can evolve with a spatial 
structure that develops as a scale-invariant system exhibiting fractal self-similar structure. In 
statistical mechanics critical phenomena can occur with transition points. Non-equilibrium 
systems undergo phase transitions with attractors. But in these dynamical systems the 
critical point can be reached only by a parameter fine-tuning and therefore takes on an 
accidental description of nature. Inflationary models, as mentioned previously, fall into this 
catagory of physical systems. But physical systems exist which evolve as self-organized criti- 
cal structures independent of the initial conditions and with no fine-tuning of the parameters 
involved in the system. A well-known example is the sand-pile model studied by Bak, Tang 
and Wiesenfeld [§J. The sand-pile model can also be pictured as a system of coupled damped 
pendula. Energy is dissipated at all length scales. 

We shall explore our model of the universe using the Lemaitre-Tolman-Bondi |5|-|8| in- 
homogeneous, spherically symmetric solution of the GR field equations. The metric is given 
by 

ds 2 = dt 2 - R' 2 (r, t)f- 2 (r)dr 2 - R 2 (r, t)dfl 2 \ (1) 

where / is an arbitrary function of r only, and the field equations demand that R(r, t) 
satisfies 

2RR 2 + 2R{1 - f 2 ) = F{r) (2) 

with F being an arbitrary function in classical GR of class C 2 ,R = dR/dt, R' = dR/dr, 
and dQ 2 = d6 2 + sin 2 Odcf) 2 . The universe is filled with dust and there are three solutions, 
depending on whether f 2 < 1,= 1, > 1 and they correspond to elliptic (closed), parabolic 
(flat), and hyperbolic (open) cases, respectively. 

Let us define the local density parameter, Q(r, t) = p(r,t)/p c (r,t). There is a correspon- 
dence between the spatial curvature and the sign of f2 — 1 : ft — 1 < 0, f 2 > 1 (open), 



Q — 1 = 0, f 2 = 1 (flat), Q — 1 > 0, f 2 < 1 (closed). The proper distance between two 
observers is defined by 



l-'r prop 



J ' Bf{r,t)f-\r)dr, (3) 



and 



L ±piop = j R(r,t)dn. (4) 

Let us now turn to a study of our one- dimensional model of the universe. Our model 
universe has a large number of metastable states, which grows exponentially with the ex- 
pansion of the universe. These metastable states correspond to attractors associated with 
an array of n stable fixed points. The universe evolves towards a minimally stable state. Let 
us use a simple model comparable to the "sand-pile model" [§]]. The numbers fl n denote the 
coupled values of the density profile associated with the forces of expansion counter-acting 
gravity. Then we have 

Q n = E(n)-E(n + l) (5) 

between successive values of Q during the expansion. 

The Q n satisfy a nonlinear, discretized diffusion equation with a threshold condition. We 
have 

fi n -> n n + 1, (6a) 

fi n _! -► fi„_! - 1. (6b) 

When Q n > Q c , where Q c is a fixed critical value of Q, a unit of density profile randomly 
moves down: 

n n± i -> Q n ±i + 1, for Q n > Q c . (7b) 

Eventually all the fl n reach the critical value Q = Q c , which is the least stable of the 
stationary states. Any additional Q n simply falls back from site to site, falling off the end 



n = N and leaving the universe in the minimally stable state. This is the cellular automaton 
model of the universe where the state of the discrete variable Q n at time t + 1 depends 
cooperatively on the state of the variable and its neighbors at time t. The dynamics leading 
to the least stable state at some time t = t s in the expansion is completely independent 
of how the universe began. We can randomly add expansion "slope" , Q n — > Q n + 1 and 
let the universe satisfy Eqs.(|7a|) and (ffb|). This corresponds to a universe with a random 
distribution of critical expansion differences and a uniformly increasing slope of expansion. 
Starting with a highly unstable universe with Q n > Q c for all n and letting the universe 
contract would lead to the same minimally stable state. 

For our one- dimensional universe, the stable critical state is robust to local perturbations. 
The pressure of expansion will build up to the point where the self- organized critical state 
is achieved. 

In a two-dimensional model of the universe, we have Q = Q(x,y) and [|j] 

£l(x-l,y)-*n(x-l,y)-l, (8a) 

Sl(x,y-l)->n{x,y-l)-l, (8b) 

n(x,y)^n(x,y)+2, (8c) 

and 

n{x,y)-+Q(x,y)-4, (9a) 

n{x,y±l)^>Q{x,y±l) + l, (9b) 

ft(s±l,y)->fi(z±l,y) + l for fi(x,y)>fi c . (9c) 

This corresponds to a simple Ising model of the universe with the expansion dynamics 
involving next-nearest neighbor interactions. 

In two spatial dimensions the minimal stable state is unstable against small fluctuations 
of the spacetime geometry and does not describe an attractor of the dynamics. The universe 
becomes stable at the point when the minimally stable clusters achieve the level when the 



fluctuations cannot be communicated through infinite distances, and there are no charac- 
teristic length and time scales. The power spectrum for metric fluctuations S(Ag) satisfies 
a power law 

S(Ag) ~ (Ag)-P (10) 

associated with a self-similar fractal behavior. 

We are unable to predict the self-organized critical value, Q = Q c , without further 
dynamical input from GR, but we can explain with this model how the universe has to 
evolve to a stable critical value Q = Q c , which is completely independent of the initial 
conditions and without any fine tuning of parameters. The present universe can be said to 
be at the "edge of chaos" - there is only one possible stable choice for the present expanding 
universe whatever its initial configuration. Even if the universe began in a chaotic mix- 
master state ||, it would soon self-organize itself away from chaos towards the edge of chaos 
and a stable, critical state. Observationally Q < 1 today, and this is the only possible 
state it can be in according to the fractal scaling of the expansion law and the power law 
dependence of the correlation function for the metric fluctuations, Eq. ( [TO] ) . 

According to our underlying assumption the spacetime geometry fluctuates randomly at 
some length scale f2|. If we assume that these fluctuations in the metric are very intense at 
the beginning of the universe, and that they smear out the light cones locally, then for a given 
short duration of time At after the big-bang there will be communication of information 
"instantaneously" throughout the universe which will solve the "horizon problem", and 
explain the present high degree of isotropy and homogeneity of the present universe. 

Although the picture of the evolving universe presented here is not incompatible with 
inflationary models, it goes beyond these models by providing a fundamental explanation 
for the robustness of the present universe. In contrast to inflationary models, there is no fine 
tuning of parameters required, since the attractor mechanism evolves from a self-organizing 
dynamical system with a high degree of cooperative phenomena based on a scale-invariant 
fractal structure. 



Further work must be done to incorporate the diffusion transport equation for Q with 
GR theory, and thereby provide a basis for more specific dynamical cosmological predictions 
of the present state of the universe and the evolution of galaxies. 

ACKNOWLEDGMENTS 

I thank M. A. Clayton and P. Savaria for helpful discussions. This work was supported 
by the Natural Sciences and Engineering Research Council of Canada. 



REFERENCES 

[1] A. Guth, Phys. Rev. D 23, 347 (1981); A. Linde, Rep. Prog. Phys. 47, 925 (1984). 

[2] J. W. Moffat, Stochastic Gravity, University of Toronto preprint, UTPT-96-16, |gr 



qc/96l0067. 



[3] N. J. Cornish and J. J. Levin, Phys. Rev. Lett. 78, 998 (1997); fer-qc/ 960502^ 



[4] P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. A 38, 364 (1988); P. Bak and K. Chen, 
Scientific American, January 1991, p. 46. 

[5] G. Lemaitre, Mon. Not. Roy. Astron. Soc, 91, 490 (1931); Ann. Soc. Sci. Bruxelles, A 
53, 51 (1933). 

[6] R. C. Tolman, Proc. Natl. Acad. Sci. U.S.A. 20, 169 (1934). 

[7] H. Bondi, Mon. Not. R. Astron. Soc. 107, 410 (1947). 

[8] J. W. Moffat and D. C. Tatarski, Phys. Rev. D 45, 3512 (1992); Ap. J. 453, 17 (1995). 



